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Abstract 

Based on experimental observations it is known that various biological cells exhibit a persistent random walk 
during migration on flat substrates. The persistent random walk is characterized by 'stop-and-go' movements 
: unidirectional motions over distances of the order of several cell diameter are separated by localized short 
time erratic movements. Using computer simulations the reasons for this phenomena had been unveiled and 
shown to be attributed to two antagonistic nucleation processes during the polymerization of the cell's actin 
cytoskeleton : the (ordinary) spontaneous nucleation and the dendritic nucleation processes. Whereas sponta- 
neous nucleations generate actin filaments growing in different directions and hence create motions in random 
directions, dendritic nucleations provide a unidirectional growth. Since dendritic growth exhibits stochastic 
fluctuations, spontaneous nucleation may eventually compete or even dominate, which results in a reorientation 
of filament growth and hence a new direction of cell motion. The event of reorientation takes place at instants 
of vanishing polarity of the actin skeleton. 
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I. INTRODUCTION 



Most animal cell types are 'motile' and possess the capacity to move over or through a substrate, and cell 
migration is a normal occurrence in both normal physiological cases as in the case of diseases. Unlike 
the phenomena related to 'mobility', where the collision among molecules determine their random 
movements, the random motion of motile cells requires an energy-consuming mechanism. It is now 
widely accepted that the basic engine for gliding or crawling locomotion of many living cells is the 
ATP-supported polymerization of the actin cytoskeleton. Rapidly moving cells can often move over the 
substrate without microtubules (keratocytes, neutrophils), but crawling motility always requires actin. 
One of the crucial factors for cell movement is actin polarity jl|, which correlates with a persistent ran- 
dom walk during migration. The persistent random walk is characterized by 'stop-and-go' movements 
: unidirectional motions over distances of the order of several cell diameters are separated by localized 
short time erratic movements. The main focus of the present study is the cause behind the persistent 
random walk of a cell, and hence the cause behind the correlated spontaneous change of polarity. 
It has been seen that cell motility and chemotactic migration are related to specific rearrangements of 
the actin cytoskeleton. Actin polymerizes into new filaments in regions of the cell that are specified by 
activated signalling at the plasma membrane These signals give rise to cell polarity and 

directional motility. In the absence of chemotactic stimuli, cells exhibit a persistent random walk. 
In this work we present results of our investigations on the relation between cytoskeletal actin poly- 
merization and the persistence random walk of a cell. The results are based on the analysis of Monte 
Carlo simulations of a simple model cell. 



II. MODEL AND SIMULATION TECHNIQUE 

Due to the limitation of information and computational complexity it is not possible to construct the 
model cell and it's motility including all regulatory proteins with their exact concentration for all types 
of cells. So we focus on one of the simplest cell which is the keratocyte residing in the epidermis, 
the most outer layer of the skin. The keratocyte is responsible for the formation of tissue and wound 
healing, both requiring the cell's ability to autonomously migrate within skin tissue. It is in fact one 
of the fastest moving cells |(| with a speed of 0.5 /im/s. The advantage of choosing keratocyte is that 
the crescent-like shape of the cell is almost constant as it moves 0, 0] • The motility of the keratocyte 
is assumed to be based largely on a continuously remodeling actin network. In our model we neglect 
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the cell body, i.e., the nucleus and other organelles. This cell type is known as a cytoplast 



a. 



A. The description of Membrane and Actin Models 

Membrane. The plasma membrane of a biological cell is a highly complex surface consisting of a 
lipid bilayer. The complexity of the cell membrane cannot be cap tured in our simple model membrane. 
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Il2| . the cell membrane is designed as an 



Therefore, similar to the previous successful studies 
elastic two-dimensional ring. Our model membrane is a flexible closed ring embedded on the square 
lattice. The ring is non-self-avoiding and exibit the usual random walk characteristics. Conformational 
changes of the ring are achieved by Monte Carlo methods, where a randomly chosen pair of two 
successive segments of the chain perform a kink jump or hair pin jump to one of the neighbouring sites. 
Globular Actin (G-actin). The model membrane encloses a fixed number, N, of actin molecules. 

n 

Each actin molecule has a size [6| of about 5 nm x 5 nm and is located at any of the vertices of the 
square lattice. The G-actin molecules diffuse freely from one lattice point to another. No excluded 
volume condition is imposed among G-actin molecules. But it is imposed between membrane and 
G-actins. So the membrane is impermeable for all actin molecules. 

Filament Actin (F-actin). The G-actin monomers form a rigid filament by associating with each 
other. Within a filament the actin monomers are called F-actin. According to experimental results 
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14l | . the filaments are assumed to be chemically coupled to the underlying substrate via membrane 
proteins. This attachment to the substrate provides the necessary traction forces for the advancement 
of the cell. Hence, in our model the actin filaments are assumed to be immobile. Excluded volume 
effects between filaments and diffusing G-actin molecule are neglected, whereas excluded volume effects 
among filaments and between filaments and membrane are included. 

Actin-associated Proteins. A group of actin associated proteins help the remodelling of the actin 
network, and hence control membrane protrusion. Proteins of the WASP/Scar families activate Arp2/3 
protein complexes. These complexes nucleate new actin filaments at the sides of existing filaments. 
The nascent filaments in the network arise in the form of branches from preexisting filaments. The 
newly formed filaments elongate at their ends, pushing the membrane at the leading edge forward 
until they are capped by specific proteins. For example, ADF/Cofilin creates free ends by severing 
preexisting filaments and promoting depolymerization of free filaments at their opposite ends, and 
CapZ and gelsolin cap the fast growing ends. 

Based on experimental observations it is known that Arp2/3 is activated only close to the membrane. 
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Therefore Arp2/3- induced branching processes must be expected to happen only near the membrane. 
The Arp2/3 molecule is not explicitly taken into account in our model because of its large physiological 
concentration, in particular near the cell membrane where it becomes activated. In our model, the 
width of the range of activation is taken to be 10 lattice sites. Accordingly, in our model we perfom 
branching from an existing filament only if the filament extends to the range of activation near the 
membrane. 

A typical snapshot of the simulated cell is depicted in Fig.l. 



B. Probabilities and Reaction Rates 



Polymerization. During the polymerization step the association of G-actin molecules to an existing F- 
actin filament occurs at certain rates at both ends. We adopt the same reaction rates for polymerization 



and depolymerization as used in a previous work 



121, and are listed in Table I. Due to the structural 



rotational symmetry of a single G-actin molecule, both ends of the F-actin filament are distinguishable 
: the 'barbed' (or synonymously, 'plus' or 'fast growing') end, and the 'pointed' (or synonymously, 
'minus' or 'slow growing') end. The association rates of ATP-bound G-actin to both ends of a filament 
are different , and we denote the rate constants and k£ for 'barbed' and 'pointed' 

ends, respectively. The same holds for the depolymerization rates k^ and k~ , albeit their difference is 
smaller. Since both ends of a filament are distinguishable, each F-actin filament possesses an intrinsic 
polarity. 

Nucleation. We consider two types of nucleation for F-actin. When nucleation is formed by two 
G-actin molecules, it is called 'spontaneous' nucleation which occurs with probability W n . When the 
'branching' nucleation takes place then a new filament is formed as a branch from the side of an existing 
filament. The branching nucleation is modeled following the 'dendritic nucleaction' model based on 
e„n tal o bs e™ tlonS H 13 Q. It is k now„ that the ac tl va ted proteia co mpl ex A rp2/ 3 can 
associate with an existing filament and can nucleate there a new filament as a branch from the mother 
filament at an angle of about 70°. This leads to the formation of a branched network. In our model on 
the square lattice, the branching process is implemented as follows. If a G-actin is found to be on the 
adjacent row to an existing filament, a new daughter filament is created at this site with probability 
W^ r . If there is already a filament in that row, then the nucleation attempt is rejected. Based on the 



experimental facts 



17. 



18l 3], the tip of the daughter filament is a plus end, and hence its direction of 
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FIG. 1: Two snapshot (a and b) of the model cell. The open circles represent G-actin, the filled triangles 
represent F-actin molecules. The length of the membrane is L=200, the total number of actin molecules is 
N=500. 5 



Quantity 


Value 


lattice constant a 


5 nm 


typical cell size 


100 x 100 lattice 


Monte Carlo step r 


0.875 /is 


W n 


0.01 


W br 


0.1 


W br 
or 


0.6 


w+ 

o 


1 


w b - 


0.0012 


w+ 


0.11 


w- 


0.02 



TABLE I: Model parameters and reaction probabilities. 



growth must be the same as the plus end of the mother filament. This implies, that the polarity of a 
filament, i.e., the vector connecting plus and minus ends, is determined by the spontaneous nucleation 
process and cannot be changed by dendritic nucleation. In our model, for the sake of simplicity, the 
spontaneous nucleation creates new filaments only parallel to the Y-axis, pointing at random in one of 
the two directions. 



Reaction Rules. Regarding the reaction rules we follow our earlier work [12 1. In brief, the reaction 
rules are as follows. First we randomly choose one actin molecule. Then for the selected actin molecule, 
a random choice with a equal probability is made between an association or dissociation process. Then 
four cases may happen. 1) If for a G-actin molecule the dissiocation process is selected, the step is 
stopped. 2) If the accociation process is selected for G-actin molecules, then reaction may happen 
provided the neighbour G or F-actin exist and the probability W n or W£. or W£ or fullfill the 
condition W > rj, where rj is a random number < r\ < 1. 3) If the association process for an F-actin 
is chosen, then a successful reaction may happen provided a G-actin as a neighbor exists. 4) If a 
dissociation process of a F-actin molecule is selected, then the process will take place with probability 
W b - or W~ or W£. . 

Reaction Probabilities and Monte Carlo steps. Since the experimental reaction rates k bp cannot 
be used directly in simulations, but instead reaction probabilities Wjj 2 , one must establish a relation 
between them. In order to calculate these probabilities, we define W b + = 1 for barbed end polymerization 
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and calculate W£ for the pointed end using the relation 

W+ k 



-t = -t, (2-1] 
W+ k+ ' 



The same can be done for the depolymerization rates W^~ p , which we choose W~=0.02 and calculate 
by using the corresponding relation to Eq. (|2.1|) . The value of has been choosen in order to 
ensure that the filaments are much shorter than the cell size during the simulations. A basic question is 
how many Monte Carlo steps, m, for Brownian motion membrane and G-actins have to be performed 
between two sucessive chemical rection attempts. Based on the experimentally known reaction rates 



and the diffusion coefficient of G-actins, a detailed analysis 
value. All model parameters are summarized in Table I. 

III. RESULTS AND DISCUSSIONS 
A. Persistent Random walk 



12j have shown that m = 5 is an appropriate 
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it is known that Arp2/3- 



From experimental observations 0, 0, |3| and from simulation studies 
induced dendritic nucleation is essential for the appearance of a persistence random walk (PRW) of 
certain cells. The PRW is characterized by 'stop-and-go' movements : unidirectional motions with 



almost constant velocity over distances of the order of several cell diameters are separated by localized 
short time erratic movements. After each interruption of the ballistic motion, the cell continues to 
move in a different direction. A typical trajectory of the center of our model cell is shown in FigJ2] as 
an example from our simulations. The persistency of the cell's random walk can be deduced from the 
time-dependent mean square displacement of the cell 

R 2 (t) = ([Y(t)-Y(0)] 2 ) (3.1) 

which is shown in FigJSJ The dynamics at short and long times are governed by diffusion, whereas for 
intermediate times, the cell motion exhibits a unidirectional drift during a certain persistence time. The 
unidirectional motion has been explained as a consequence of the autocatalytic nature jsiL 22^ of the 
dendritic polymerization kinetics. The branching processes induce an explosive growth of filaments in 
a certain direction which leads to membrane protrusions and concomitant cell movements in the same 
direction. This motion lasts only for a certain persistence time, then the direction of motion changes. 
In the following subsection we discuss the reasons for this dynamic instability. It should be noted that 
insights to the physical basis of the cell's persistency would support, in particular, our understanding 
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FIG. 2: The trajectory Y(t) (full line, left scale) and polarity P(t) (dotted line, right scale) of a cell exhibiting 
a persistent random walk, as function of time t in units of Monte Carlo steps, MCS. 

of the physical basis of chemotaxis, i.e., the prolongation of the persistency of the cell's motion by an 
external signal. 

B. Dynamic Instability 

The fundamental quantity which governs the migration of the cell is its polarity and its mechanism of 
self-perpetuation. The cell polarity is macroscopically and dynamically characterized by the formation 
of leading and trailing membrane edges which are a consequence of the cell's coupling to the polarity 
of the enclosed actin cytoskeleton. The polarity of the cytoskeleton itself is the average orientation 
of its filaments, where each actin filament possesses a polarity due to its plus and minus ends. This 
spontaneous internal polarity of the actin network determines the direction of cell motion and is main- 
tained without external signals (e.g. chemotaxis) only for a certain time (persistent random walk). 
Specifically, the formation and regulation of cell polarity is achieved by a complex protein signaling 
network with positive and negative feedback loops. We do not discuss the details here which are beyond 
the scope of this paper. We note that different signals converge, among others, on the activation of an 
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FIG. 3: Typical mean square displacement of a cell exhibiting a persistent random walk. 



Arp2/3 protein complex, which leads to branching and the autocatalytic polymerization. 

In our model, we define the polarity of the cell by the polarity of the total actin network, i.e., by the 

difference in the number of filaments pointing in opposite directions, 

p _ {F up ~ Fdn) ^ 2^ 

Ftotal 

where F up is the number of filaments pointing in positive Y-direction, Fd n in negative Y-direction, and 
Ftotai is the total number of filaments. The range of the polarity parameter is — 1 < P < 1. The 
trajectory of this quantity is shown in Fig. El (dotted line, right scale). From this result one observes 
that at certain times, where P(t) = 0, the cell changes its direction of motion to the opposite as 
indicated by the corresponding displacements Y(t). 

After the analysis of several simulations under different conditions of the model cell, it must be concluded 
that the cause behind this phenomena is the existence of two antagonistic processes, the spontaneous 
and the dendritic nucleation processes. 

During a spontaneous nucleation process, where two G-actin monomers form a new F-actin filament, 
the orientation of the new filament is determined at random. This is in contrast to the dendritic 
nucleation process, where the daughter filament branches off the mother filament and therefore always 
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adopts the same growth direction as the mother filament. Hence, spontaneous nucleation changes the 
sign of polarity, whereas dendritic nucleation increases the absolute value of polarity. The first leads 
to random motion, whereas the latter prolongs unidirectional motion. Since the polarity fluctuates in 
time, the spontaneous nucleation may eventually dominate and causes a change of polarity, and hence 
a change of vectorial motion. 

This view is corroborated by simulations of various limiting cases, which are shown in Fig0] 
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FIG. 4: Mean square displacements of a cell under different conditions and at different probablilities of sponta- 
neous nucleation, wn. 's' represents spontaneous nucleation, 's+d' represents both spontaneous and dendritic 
nucleation, and 'd' represents dendritic nucleation only. 

In one limiting case, the cell motion is caused solely by spontaneous nucleations. The corresponding 
curve is shown in FigJUby the full lines with label 's'. At two different nucleation probabilities, W n , 
the mean square displacements indicate ordinary diffusion, R 2 (t)/t ~ const., at almost all times. This 
is expected due to the random orientation of newly established filaments. A very weak persistency may 
be expected at small nucleation rates. 

In the other limiting case, the cell motion is caused solely by dendritic nucleation. Since in this case 
the orientation of the filaments is preserved, the corresponding curve (dotted line with label 'd') has 
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to exhibit a clear drift, R 2 (t)/t ~ t, which takes place for times larger than a characteristic time Td ~ 
7 x 10 3 . The characteristic time scale Td separates periodes of cell advancements and localized cell 
displacements ('resting'). Since the activation range of Arp2/3 extends only a few lattice sites from the 
membrane, the rapid advancement of the leading edge of the membrane together with the rapid auto- 
catalytic processes of branching near the leading edge deprive this area of available G-actin molecules 
which are necessary for continuation of these processes. Therefore, actin network and cell remain tem- 
porarily, during t < Td, essentially at rest. While recruiting G-actins by diffusion from the minus to 
the plus ends of the filaments near the activation zone, the membrane of the cell performes random 
displacements around the immobile filament network, which leads to time-independent displacements 
of the cell, R 2 (t)/t ~ l/t at t < Td. It is still unknown whether this result from simulations is related 
to the experimentally observed |23( excitation waves of F-actin assembly near the membrane and the 
correlated cell advancements. 

The interesting intermediate case, when a few events of spontaneous nucleation interfere with den- 
dritic nucleations, is shown in Fig0] by the broken line with label 'd+s'. At a nucleation ratio of 
W n /W£. =0.015, again a drift of the cell is observed albeit only up to a typical persistence time 
t p m 10 6 . At larger times, ordinary diffusion R 2 (t)/t = const, is observed indicating the dominance of 
spontaneous nucleation processes. For t < t p , the mean square displacements are qualitative the same 
as for the dendritic case 'd'. 

IV. SUMMARY AND CONCLUSIONS 

The dynamics of a motile cell with unidirectional and bidirectional branching is studied with tread- 
milling actin polymerization kinetics. The filaments growing on both sides with treadmilling rates and 
branching in both directions leads to a persistent random walk whose persistency becomes infinite and 
leads to unidirectional drift in the case of one side branching. The dynamics are governed by diffusion 
at short times followed by drift at longer times. 

The phenomena is shown to be attributed to two antagonistic nucleation processes during the poly- 
merization of the cell's actin cytoskeleton : the (ordinary) spontaneous nucleation and the dendritic 
nucleation processes. Whereas spontaneous nucleations generate actin filaments growing in different 
directions and hence create motions in random directions, dendritic nucleations provides a unidirec- 
tional growth. Since dendritic growth exhibits stochastic fluctuations, spontaneous nucleations may 
eventually compete or even dominate, which results in a reorientation of filament growth and hence a 
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new direction of cell motion. The event of reorientation takes place at instants of vanishing polarity of 
the actin skeleton. 
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